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Resummed jet rates with heavy quarks in e+e colUsions * 
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Theory Division, CERN, CH-121 1 Geneva 23, Switzerland 

Expressions for Sudakov form factors for heavy quarks are presented. They are used to construct resummed jet rates in e^e~ annihila- 
tion. Predictions are given for production of bottom quarks at LEP and top quarks at the Linear Collider. 



1 Introduction 

The formation of jets is the most prominent feature of per- 
turbative QCD in e+e^ annihilation into hadrons. Jets 
can be visualized as large portions of hadronic energy or, 
equivalently, as a set of hadrons confined to an angular 
region in the detector. In the past, this qualitative defini- 
tion was replaced by quantitatively precise schemes to de- 
fine and measure jets, such as the cone algorithms of the 
Weinberg-Sterman [ 1 j type or clustering algorithms, e.g. 
the Jade [ 2 ^ or the Durham scheme (fc^ scheme) [ 4 1. A 
refinement of the latter one is provided by the Cambridge 
algorithm ["S^. Equipped with a precise jet definition the 
determination of jet production cross sections and their in- 
trinsic properties is one of the traditional tools to investi- 
gate the structure of the strong interaction and to deduce 
its fundamental parameters. In the past decade, precision 
measurements, especially in e+e^ annihilation, have es- 
tablished both the gauge group structure underlying QCD 
and the running of its coupling constant over a wide 
range of scales. In a similar way, also the quark masses 
should vary with the scale. 

A typical strategy to determine the mass of, say, the 
bottom-quark at the centre-of-mass (cm.) energy of the 
collider is to compare the ratio of three-jet production cross 
sections for heavy and light quarks [ '6"7"8"9"r0 1 . At jet 
resolution scales below the mass of the quark, i.e. for glu- 
ons emitted by the quark with a relative transverse momen- 
tum k± smaller than the mass, the collinear divergences 
are regularized by the quark mass. In this region mass 
effects are enhanced by large logarithms ln(mb/fcx), in- 
creasing the significance of the measurement. Indeed, this 
leads to a multiscale problem since in this kinematical re- 
gion also large logarithms ln(^/s/fcj^) appear such that 
both logarithms need to be resummed simultaneously. A 
solution to a somewhat similar two-scale problem, namely 
for the average sub-jet multiplicities in two- and three-jet 
events in e+e" annihilation was given in [ 11|. We re- 
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port here on the resummation of such logarithms in the fcj^- 
like jet algorithms [ 12| and provide some predictions for 
heavy quark production. A preliminary comparison with 
next-to-leading order calculations of the three-jet rate [ 
|6|[l3l[ll[i5J is presented. 



2 Jet rates for heavy quarks 

A clustering according to the relative transverse momenta 
has a number of properties that minimize the effect of 
hadronization corrections and allow an exponentiation of 
leading (LL) and next-to-leading logarithms (NLL) [ "41 
stemming from soft and collinear emission of secondary 
partons. Jet rates in k±_ algorithms can be expressed, up 
to NLL accuracy, via integrated splitting functions and Su- 
dakov form factors [ 4|. For a better description of the jet 
properties, however, the matching with fixed order calcu- 
lations is mandatory. Such a matching procedure was first 
defined for event shapes in [ 16J. Later applications include 
the matching of fixed-order and resummed expressions for 
the four-jet rate in e+e" annihilation into massless quarks [ 
Q? 18 1. A similar scheme for the matching of tree-level 
matrix elements with resummed expressions in the frame- 
work of Monte Carlo event generators for e'^e^ processes 
was suggested in [ .19J and extended to general colhsion 
types in r l20l . 



We shall recall here the results obtained in ['12| for heavy 
quark production in e+e" annihilation. In the quasi- 
collinear limit [ 12111221 . the squared amplitude at tree-level 
fulfils a factorization formula, where the splitting functions 
Pab for the branching processes a ^ 6 + c, with at least 
one of the partons being a heavy quark, are given by 



Pqq{z,<i) = Cf 
PgQiz,q) = Tr 



l + z^ 2z(l -z)m2 



1 — z g2 + (1 — z)2m^ 

„ N 2z(l — z)m 

l-2z{l-z) + ' ^ ' 



(1) 



where z is the usual energy fraction of the branching, and 
is the space-like transverse momentum. As expected. 
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these splitting functions match the massless spUtting func- 
tions in the hmit to — > for fixed. The splitting function 



z 1 — z 

+ + z{l-z) 

1 — z z 



(2) 



obviously does not get mass corrections at the lowest order. 
Branching probabilities are defined through [ 1121 
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with 



rQ(g,g,m = 0) = 2Cp(ln^-^ 



(4) 



and the Sudakov form factors, which yield the probabiUty 
for a parton experiencing no emission of a secondary par- 
ton between transverse momentum scales Q down to Qq, 
read 



Aq(Q,(3o) = exp 
A<,(Q,Qo) = 
exp 



Q 



dq asiq) 
q TT 
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Qo 
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q TT 
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AfiQ; Qo) = [Aq(Q, Qo)] /Ag(Q, Qo) , 



(5) 



where rf{Q,q) accounts for the number n^''''' of active 
light or heavy quarks. Jet rates in the k± schemes can be 
expressed by the former branching probabilities and Su- 
dakov form factors. For the two-, three- and four-jet rates 

7^2 - [Aq(q,Qo)]', 
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where Q is the cm. energy of the colliding e+e^, and 
Qo = J/cutQ^ plays the role of the jet resolution scale. 
Single-flavour jet rates in Eq. (|6} are defined from the 
flavour of the primary vertex, i.e. events with gluon split- 
ting into heavy quarks where the gluon has been emitted off 
primary light quarks are not included in the heavy jet rates 
but would be considered in the jet rates for light quarks. 

In order to catch which kind of logarithmic corrections are 
resummed with these expressions it is illustrative to study 
the above formulae in the kinematical regime such that 
Q ^ TO ^ Qo. Expanding in powers of as, jet rates 
can formally be expressed as 



E 

fc=n-2 



2^^ 



kl ' 



(7) 



(n) 

where the coefficients cj.^ are polynomials of order I in 
Ly = ln(l/2/cut) and L,„ = ln(TO2/QQ). The coefficients 
for the first order in as are given by 



-12 



„(2) 
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For second order q^, with n active flavours at the high 
scale, the LL and NLL coefficients read 



J2) 
^24 
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An approximate way of including mass effects in massless 
calculations, that is sometimes used, is the "dead cone" [ 
|23J approximation. The dead cone relies on the observa- 
tion that, at leading logarithmic order, there is no radiation 
of soft and collinear gluons off heavy quarks. This effect 
can be easily understood from the splitting function Pqq 
in Eq. Q. For g <C (1 — z) m this splitting function is not 
any more enhanced at z 1. This can be expressed via 
the modified integrated sphtting function 
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Terms ^ (/?„ — /3„_i) in the NLL coefficients, where the 
/3-function /3„ for n active quarks is given by 



Pn = 



lie A- 'in 
12 



(10) 



are due to the combined effect of the gluon splitting into 
massive quarks and of the running of as below the thresh- 
old of the heavy quarks, with a corresponding change in the 
number of active flavours. With our definition of jet rates 
with primary quarks the jet rates add up to one at NLL ac- 
curacy. This statement is obviously realized in the result 
above order by order in as- 

The corresponding massless result [ "41 is obtained from 
Eqs. (|8} and (|9} by setting L„i 0. Notice that Eqs. (|8j 
and (|9jl are valid only for m ^ Qq and therefore m — > 
does not reproduce the correct limit, which has to be 
smooth as given by Eq.(|6jl. Let us also mention that for 
Q > TO there is a strong cancellation of leading logarithms 
and therefore subleading effects become more pronounced. 



rQ(Q,g,TO = 0) 
2Cp\n(^) Q(n 



9) 



(11) 



To obtain this result the massless splitting function has 
been used, which is integrated with the additional con- 
straint z > 1 — q/m. We also compare our results with 
this approximation. 



3 Numerical results and comparison with 
fixed order calculations 

The impact of mass effects can be highlighted by two ex- 
amples, namely by the effect of the bottom quark mass in 
e"'"e~ annihilation at the Z-pole, and by the effect of the 
top quark mass at a potential Linear Collider operating in 
the TeV region. 

With TOf, = 5 GeV, Mz = 9L2 GeV, and = 
0.118, the effect of the &-mass at the Z-pole on the two- 
and three-jet rates is depicted in Fig. ^ (left). The result 
obtained in the dead cone approximation is shown in Fig. [J 
(right). Clearly, by using the full massive sphtting func- 
tion, the onset of mass effects in the jet rates is not abrupt 
as in the dead cone case and becomes visible much ear- 
lier. Already at the rather modest value of the jet resolu- 
tion parameters of j/cut — 0.004, the two-jet rate, includ- 
ing mass effects, is enhanced by roughly 4% with respect 
to the massless case, whereas the three-jet rate is decreased 
by roughly 3.5%. For even smaller jet resolution param- 
eters, the two-jet rate experiences an increasing enhance- 
ment, whereas the massive three-jet rate starts being larger 
than the massless one at values of the jet resolution param- 
eters of the order of ?/cut ~ 0.001. The curves have been 
obtained by numerical integration of Eq. (|6|l. Furthermore, 
in order to obtain physical result the branching probabili- 
ties have been set to one whenever they exceed one or to 
zero whenever they become negative. 

While in the case of bottom quarks at LEPl energies the 
overall effect of the quark mass is at the few-per-cent level, 
this effect becomes tremendous for top quarks at the Linear 
ColHder (Fig.|3. 

In Fig. |3l leading order (LO) and next-to-leading order 
(NLO) predictions for three-jet rates are compared with 
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Figure 1. Effect of a b-mass of 5 GeV on the single-flavour two- and tiiree-jet rate at LEPl energies as a function of the jet resolution 
parameter in the fcx scheme. In the left plot this effect is treated through the full inclusion of masses into the splitting function, see 
Eq. J3j' whereas in the plot on the right hand side this effect is modeled through the dead cone, see Eq. GI}. 
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Figure 2. The effect of a t-mass of 175 GeV on the single-flavour two- and three-jet rate as a function of the jet resolution parameter, 
at a potential Linear Collider operating at cm. energies of 1 TeV. Again, in the left plot this effect is treated through the full inclusion 
of masses into the splitting function, see Eq. whereas in the plot on the right hand side this effect is modeled through the dead cone, 
see Eq. <1 1> . 
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Figure 3. Comparison of LO, NLO and NLL predictions for the single-flavour three-jet rate as a function of the jet resolution parameter 
in the Cambridge algorithm for bottom quark production at LEPI energies (left plot), and top quark production at a potential Linear 
Collider operating at cm. energies of 1 TeV (right plot). 



the NLL result showed in the previous plots. Fixed or- 
der predictions for 6-quark production clearly fail at very 
low values of j/cut, by giving unphysical values for the jet 
rate, while the NLL predictions keep physical and reveal 
the correct shape. The latter is an indication of the neces- 
sity for performing such kind of resummations. Fixed or- 
der predictions work well for top production at the Linear 
Collider, a consequence of the strong cancellation of lead- 
ing logarithmic corrections, and are fully compatible with 
our NLL result. 



4 Conclusions 

Sudakov form factors involving heavy quarks have been 
employed to estimate the size of mass effects in jet rates in 
e^e~ annihilation into hadrons. These effects are sizeable 
and therefore observable in the experimentally relevant re- 
gion. A preliminary comparison with fixed order results 
have been presented, and showed good agreement. Match- 
ing between fixed-order calculations and resummed results 
is in progress [|24|. 
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